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Abstract -Synchronization problems in complex networks are very often studied by researchers 
due to its many applications to various fields such as neurobiology, e-commerce and completion of 
tasks. In particular, Scale Free networks with degree distribution P{k) ^ k~^, are widely used in 
research since they are ubiquitous in nature and other real systems. In this paper we focus on the 
surface relaxation growth model in Scale Free networks with 2.5 < A < 3, and study the scaling 
behavior of the fluctuations, in the steady state, with the system size N. We find a novel behavior 
of the fluctuations characterized by a crossover between two regimes at a value of A = N* that 
depends on A: a logarithmic regime, found in previous research, and a constant regime. We propose 
a function that describes this crossover, which is in very good agreement with the simulations. We 
also find that, for a system size above N*, the fluctuations decrease with A, which means that the 
synchronization of the system improves as A increases. We explain this crossover analyzing the 
role of the network’s heterogeneity produced by the system size N and the exponent of the degree 
distribution. 


Since a great variety of systems can be represented by complex networks, over the last 
decades many researchers have studied both the topology and processes that evolve on top 
of these networks. Systems such as neural networks, the Internet and airlines networks m 
can be described by a set of nodes connected by links that represent a relationship between 
them, such as an electric impulse, friendship or air traffic. Many of these real networks were 
found to be characterized by a Scale Free (SF) topology, given by a degree distribution 

P{k)^k-\ ( 1 ) 

where k is the degree of the nodes and m < k < kmax , where m and k^ax are the minimum 
and maximum degree respectively, and A represents the broadness of the distribution. On 
most real systems, such as the World Wide Web or metabolic networks, it was found that 
2 < A < 3 |T|[^. 

More recently, research has focused on dynamical processes taking place on the under¬ 
lying network [4lII0). Particularly, many mathematical and numerical models have been 
elaborated to study the problem of synchronization [nHn, a phenomenon present in the 
behavior of many collective systems. In these processes the state of the system evolves to a 
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synchronized state, where the coupled units adjust their dynamics with one another. Exam¬ 
ples of synchronization can be seen in brain processes HZ] or data distribution [TBfED] : in a 
network made up of processors that distribute the task load, the system is best synchronized 
when the process minimizes the waiting time in each processor. For these kinds of systems, 
a scalar field h is usually defined on the network and it is of interest to measure the fluctu¬ 
ations of h. This problem can be studied mapping it into a non-equilibrium surface growth 
problem, where the scalar field /li, with i = 1, ...N and N is the system size, represents the 
“height” of the node i, and the fluctuations, also called roughness of the system, are given 
by 


and 


W{t) 




1 ^ 


1 

(H*)) = 


( 2 ) 

(3) 


is the mean value of h at time t. The roughness has two regimes, one where W{t) increases 
until it saturates at a constant value Ws in the second regime, which depends on the topology 
of the system. In Euclidean lattices one of the most studied equations on surface growth 
is the Edward Wilkinson (EW) equation [21], which belongs to the same universality class 
that the stochastic growth model of surface relaxation to the minimum (SRM) [51]. In this 
model, a node of the network is randomly selected and a “particle” is placed on the node 
with the lowest height among the selected node and its neighbors. According to this rule, 
the nodes with higher heights distribute their excess of particles to their neighbors with 
lower heights. Pastore y Piontti et al. [12] studied the SRM model on SF networks through 
numerical simulations and found that, on the steady state, the behavior of the fluctuations 
with the system size N is given by 


J const. , for A > 3 ; 
In N , for A < 3 . 


(4) 


Unlike Euclidean lattices, in complex networks one cannot extend the discrete nature of the 
network to a continuum, and thus the dynamics of the network are not well represented by 
a continuum equation such as the EW equation. However, using a discrete Laplacian and a 
mean field approximation, Korniss et al. [13] and Guclu et al. [23] found that in the limit 
N ^ oo, Ws increases with A. Solving the discrete EW equation numerically for finite size 
systems, in [12] the authors found that Wg decreases with N, which is not representative of 
any growth model. With a different approach, La Rocca et al. m developed a Langevin 
stochastic equation that describes the evolution of the interface, and solved it up to second 
order by numerical integration for finite system sizes, recovering Eq. O- 

In this letter we mainly consider SF networks with A < 3 because they are representative 
of abundant systems in nature. We find that, for the SRM model, Ws has a crossover from 
a logarithmic regime to constant regime, at a characteristic value of N that depends on the 
topology of the network. Also, we find that, for system sizes above this characteristic size, 
Ws increases as A decreases. 


By stochastic numerical simulations, we study the SRM process on SF networks. To 
generate the network, we use the Molloy Reed algorithm (MR) or configurational model 
|24j and we use a minimum degree m = 2 because for this value we have a high probability 
of obtaining only one component and thus a single interface [25] . As for the maximum 
degree, the network has a natural cut-off given by ^ [53] and no structural cut-off 
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is imposed, although this alternative will be discussed later. We define a scalar field h on 
the network, which represents the system’s feature we want to study, so that each node is 
assigned a height hi, with i = 1,N. At t = 0, we allocate the nodes with a random height 
between 0 and 1. This initial condition does not affect the scaling behavior of the roughness 
in its steady state. At each time step, we deposit a particle on a node randomly selected 
with a probability 1/A^. Denoting the selected node by i and the set of its ki neighbors by 
Vi, we simulate the SRM process according to the following rules 

1) if hi < hjV j e Vi ^ hi = hi + 1 , 

2) else, I G Vi : hi < hi and hi < hj\/j ^ I, j evi ^ hi = hi + I , 

and compute the roughness at the saturation as a function of the system size. 




Fig. 1: a) Ws as a function of A in a linear-log scale, with A = 2.6 (O); 2.8 (□), 3.0 (*). We 
observe two regimes so that Ws behaves as ~ InN for N smaller than a certain characteristic size 
N*, and it increases at a slower rate for N > N*, finally becoming constant for a sufficiently large 
N. b) Wa as a function of A in a linear-log scale for A = 2.6. ft can be seen how we obtain the 
characteristic system size A* for the intersection of the two regimes. 


In Fig. [l](a) we plot the fluctuations Wg as a function of N for different values of A. We 
can see that for A = 3 we obtain the behavior predicted by Eq. O, i.e, the fluctuations go 
rapidly to a constant when the system size increases. For A < 3 and for a range of system 
sizes (A < 10^), the behavior of Ws is logarithmic, which agrees with Eq. (|4]), but when 
N increases, the fluctuations increase slower than a logarithmic, reaching a constant that is 
independent of N and only depends on A. The scaling behavior of Ws suggests that above a 
certain system size N = N* (which depends only on A), the fluctuations become independent 
of N and therefore the system has the same degree of synchronization. It is worth noticing 
that the second regime was not seen in [T^] since in their research they simulated systems 
smaller than A*. 

We then estimate the system size for which the behavior of the fluctuations changes from 
a logarithmic regime to a constant, i.e., the crossover between these two regimes. In Fig. [T] 
(b) we plot Ws as a function of A only for A = 2.6 in order to show how A* is determined. 
We compute A* for different values of A, and we see that A* decreases with A for A < 3, 
and A* 0 for A > 3, as can be seen in Fig. [2l 

The behavior of Ws with A can be described as follows 

f b ln(N) , for A < A* ; 

W~ ,tor]V>W, 

where b = 6(A) and W^ = IFg°“(A) is the roughness value in the thermodynamic limit 
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Fig. 2: N* as a function of A, for SF networks with 2.6 < A < 2.9 with m — 2. 


(above TV*). Thus we propose a scaling function f{N/N*) where 
f{N/N*) ^ 

Then, the behavior of Wg for all the values of N can be expressed as 


f ln{N/N*) , for N/N* < 1; 
\ 0 , for N/N* > 1. 


Wg = + b f{N/N*) . 


To lose all dependency with A, we work with the expression (Ws°° — Wg^/h so that 


(irr - Wg)lh - -f{N/N*) 


-ln{N/N*) , for N/N* < 1 

0 , for N/N* > 1 


In Fig. [3] we plot (W^°“ — Wg)/b as a function of N/N*. From the plot we can see that 
the curves indeed overlap, which shows that our scaling hypothesis is correct. 



N/N 


Fig. 3: (Ws°° — Wa)/h as a function of N/N* in a linear-log scale for A = 2.6(o), 2.7(<l) and 2.8(0)- 
The curves overlap, which proves the scaling hypothesis given by Eq. ([S]) correct. 


To understand the dynamics of the system as N increases, we study the behavior of the 
fluctuations relative to the topology of the network, specifically the degree of the nodes. 
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We compute the mean height of the nodes with degree k, denoted by hk- In Fig. |4] we 
show hk — (h) as a function of fc for A = 2.6 and different values of N, below and above 
N*. We see that, in average, the difference between the height of the nodes and the mean 
value of the entire network, increases with k. Thus, for the SRM model, nodes with high 
degree worsen synchronization while small degree nodes improve it. This is because the 
nodes with high degree receive the excess of particles of their low degree neighbors, which 
are the majority in SF networks. In order to understand the effect of the high degree nodes 
on the behavior of Ws , we study the SRM process for a network with a structural cut-off of 
ks ~ [SSIET], which is smaller than kc for A < 3 (not shown here). For this case, we 
also find a crossover between two regimes at a characteristic size, although it is larger than 
the one found with no structural cut-off. This means that the hubs contribute to the finite 
size effects for << N* with A < 3. 

We also notice that the behavior of hk — (h) does not depend on N for small k, which 
means that low degree nodes do not contribute to finite size effects on Wg- However, for 
larger fc, hk — {h) increases with N. As k increases, the rate of increase of hk — (h) decreases, 
so that the nodes have heights more similar to one another. This behavior combined with 
the fact that the probability of high connectivities is very low determines that nodes with 
high degree do not contribute to an increase in the fluctuations for N » N*. Due to this 
combined effect Wg approaches to a constant for N >> N*. On the other hand, as A 
increases, the ratio of small degree nodes to hubs increases. This explains why, around the 
value of N*, the synchronization improves as A increases. This can be seen from Fig. [l] 
where N* decreases as A —>■ 3 and, consequently, the constant value to which Wg approaches 
gets smaller, which means that the system synchronizes better for increasing A. 



Fig. 4: hk — (h) as a function of k for N = 6044 (Q), 32768 (□), 500000 (o) and 1000000 (A), and 
A = 2.6. The inset is an enlargement for small values of fc; it can be seen that hk — {h) does not 
depend on N in this region. 


To summarize, we studied the behavior of the fluctuations Wg with the system size 
N in the steady state for the SRM model in SF networks with 2.5 < A < 3, aiming to 
understand the role of finite size effects in the system’s synchronization. We found a crossover 
between two different regimes at a characteristic size N* below which Wg has a logarithmic 
dependence with N and above which Wg is constant. We measured N* for different values 
of A and found that N* decreases with A. We also found that the synchronization enhances 
as A increases. The behavior of Wg with N and A is determined by the performance of high 
degree nodes in the dynamics of the system, which reach heights above the mean value and 
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worsen the synchronization. However, the rate of increase of the heights compared with the 
average height decreases with k and, given that the degree distribution also decreases with fc, 
we conclude that high degree nodes do not contribute to an increase of Wg for N ^ N*. It 
is important to mention that eventhough high degree nodes are the responsible of the finite 
size effects observed for N < N*, the explanation of the logarithmic behavior of Wg goes 
beyond the aim of our actual research. However, this will be the scope of future researches. 
As for the heterogeneity of the network, as A increases, the proportion of high degree nodes 
decreases and the previous effect is noted for smaller system sizes. In the limit A —>■ 3, 
N* —>■ 0 and Wg is constant for all N. 
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